Quantum critical transition from charge-ordered to superconducting state in the 
triangular lattice negative-U extended Hubbard model 
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We demonstrate a robust frustration-driven charge-order to superconductivity transition in the 
half-filled negative-U extended Hubbard model. Superconductivity extends over a broad region of 
the parameter space. We argue that the model provides the correct insight to understanding uncon- 
ventional superconductivity in the organic charge-transfer solids and other quarter-filled systems. 
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Spatial broken symmetries such as antiferromagnetism 
(AFM) and charge ordering (CO) are proximate to 
superconductivity (SC) in a number of exotic sys- 
tems, including the cupratesi, Na a: Co0 2 ■ yH 2 02, /3- 
Nao.33V 2 Og£ and organic charge-transfer solids (CTS) 
such as (TMTCF) 2 X (here C=S or Se, and X are closed- 
shell anions)^, (BEDT-TTF) 2 X (hereafter ET 2 X)i and 
EtMe 3 Z[Pd(dmit) 2 ] 2 , Z = P, As 5 . Unlike in the cuprates, 
SC in the CTS is reached not by doping, but on appli- 
cation of hydrostatic or uniaxial pressure. CTS crys- 
tals often consist of anisotropic triangular lattices of 
dimers of the active cationic or anionic molecules. The 
average number of charge carriers n per molecule is 
indicating that n per dimer unit cell is 1. Since 
1 triangular lattice provides the classic tem- 
plate for the resonating valence bond (RVB) electronic 
structure within the nearest-neighbor (n.n.) Heisenberg 
Hamiltonian^, the idea that spin frustration drives an 
AFM-to-SC or spin liquid-to-SC transition in the CTS 
has acquired popularity 7 . Within this picture, pressure 
makes the effective anisotropic n = 1 triangular lattice 
more isotropic, and SC occurs over a narrow range of 
anisotropy between the more robust AFM and the para- 
magnetic metal (PM). Numerical quantum Monte Carle ^ 
and exact diagonalization- calculations, however, have 
failed to find superconducting correlations in the the tri- 
angular lattice repulsive n = 1 Hubbard model, casting 
doubt on the mean-field techniques that find SC within 
the Hamiltonian. Experimentally, the situation is com- 
plex, (i) SC in certain CTS is proximate to CO instead 
of AFMi£. This has led to yet other mean-field models 
with additional Coulomb parameter of charge-fluctuation 
mediated SGAi. (ii) The insulating phase proximate to 
SC in EtMe 3 [Pd(dmit) 2 ] 2 is not an AFM but a valence 
bond solid (VBS), with charge disproportionation be- 
tween molecules^, (iii) AFM is missing in the insulating 
state of k— ET 2 Cu 2 (CN)3 with a nearly isotropic trian- 
gular lattice^. There occur inhomogeneous charge local- 
ization and sharp decrease in spin susceptibility below 10 
K-ii, which may also be signatures of a static or fluctuat- 



ing VBS-like state. Whether or not frustration can drive 
transition to SC from an ordered state therefore remains 
an open and intriguing question. 

In the present paper we demonstrate a robust 
frustration-driven SC within the n = 1 negative-?/ ex- 
tended Hubbard model (EHM) with n.n. Coulomb re- 
pulsion V. Although the literature on the negative-?/ 
Hubbard Hamiltonian is vast, the model has been inves- 
tigated primarily for bipartite lattices. With repulsive 
V, there is strong tendency to CO in bipartite lattices 
and SC is absentAS. Frustrated lattices have been inves- 
tigated within the Hamiltonian for V = onlyi&. Now 
SC dominates and CO is absent. Here we begin with 
the square lattice with n.n. V and electron hopping, 
when the ground state is a checkerboard CO with alter- 
nate double occupancies and vacancies in the square lat- 
tice. As the Coulomb interaction and electron hopping 
along one diagonal of the square lattice increase from 
zero, charge frustration in the emergent triangular lattice 
leads to first-order transition to a superconducting state. 
While our primary goal is to demonstrate the frustration- 
driven CO-to-SC transition, we also point out that our 
work provides insight for understanding unconventional 
SC in n = i correlated electron systems including the 
CTS. 

We consider the two-dimensional (2D) Hamiltonian, 
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on an anisotropic triangular lattice. Here c\ a creates an 
electron with spin a (|, J,) on site i, n,i t(7 — c\ a Ci ta , and 
m = X)<T ni > CT - ^ ^ s t ne on-site Hubbard interaction; (...) 
implies n.n. along the x- and y-axes, with hopping inte- 
gral and Coulomb repulsion t and V, respectively. Sim- 
ilarly, [• • • ] implies neighbors along the (x+y)-diagonal, 
with t' and V as the hopping and Coulomb integrals. 
We have chosen the same signs for t and t' , as the results 
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for n — 1 are independent of the relative sign of if. In 
the following we express all quantities in units of t. For 
simplicity we consider mostly V = V, although we have 
performed calculations also for other V . Nonzero V and 
if are both crucial for the CO-to-SC transition. 

We have performed exact diagonalizations on a 16- 
site periodic lattice. The quantities we calculate are 
the structure factor S(Q), the bond-order Bd along the 
(x + y)-diagonal, and the pair correlation function P(r), 
which are defined as, 



I'.i X > ' (2b) 

a 

P ( r ) = ]v E^T^d^+^-HVl) ( 2c ) 
i 

We begin with the numerical results for U = — 2 and 
V = V = 1, for which a first-order like CO-to-SC transi- 
tion occurs in our finite lattice at t' c = 0.5. In Figs. 1(a) 
and (b) we have plotted S(Q) versus Q for if = 0.45 and 
if = 0.5, respectively. The sharp peak at Q = (%, 7r) 
in Fig. 1(a) is a signature of the checkerboard CO, with 
doubly occupied and vacant sites alternating along the x 
and y-directions (see Fig. 4(b)). The amplitude of the CO 
remains practically unchanged between if = and 0.48 
(see below). The complete absence of the S(Q) peak in 
Fig. 1(b) indicates a sudden loss of the CO due to the 
mobility acquired by the electron pairs. The latter in 
turn leads to superconducting pair correlations, as seen 
in Fig. 1(c), where we have plotted P(r) against r (where 
r is in units of the lattice constant) for both t' — 0.45 
and 0.5. For r = 1.414, P(r) has two different values 
corresponding to diagonals along x + y and x — y, re- 
spectively. We have chosen to show the smaller of the 
two, corresponding to the nonbonded sites along x — y, 
in Fig. 1(c). The large increase in P(r) for all r (by more 
than a factor of 2 at n.n., and nearly an order of magni- 
tude at larger r) as t' changes from 0.45 to 0.5, as well 
as the weak dependence of P(r) on r beyond r = 1 are 
signatures of a SC ground state for if — 0.5. While SC 
is present in the model at electron densities 1, the 
CO-to-SC transition is unique to n = 1. Our calculations 
were done also for V' ^ V. The ground state continues 
to be CO for all if for V = 0, with the amplitude nearly 
independent of if (not shown). Furthermore, P(r) in this 
case is nearly for all r > 1. Thus both nonzero V 
and t' are essential for frustration. We have, however, 
determined that the CO-to-SC transition occurs here at 
nearly the same t' c for V' = t'. 

To illustrate the sharpness of the CO-to-SC transition 
in this case we have plotted in Fig. 2(a) S(n, n) and Bd 
against t' for the same U,V as in Fig. 1. The sudden 
drop in S(n, n) and the jump in Bd occur at the same 
t' c , indicating that the vanishing of the CO is due to sud- 
den increase in carrier mobility. Fig. 2(b) shows a plot of 




FIG. 1: (a) and (b), S(Q) vs. Q for U = -2, V = V = 1 
and t' = 0.45 and t' = 0.5, respectively, (c) P(r) vs. r, for 
t' = 0.45 (dashed curve) and t' = 0.5 (solid curve). 
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FIG. 2: (a) S(7T,7r) (filled circles) and Bd (unfilled circles) vs. 
t' for U = -3, V = V' = 1. (b) P(r max ) vs. t' for the same 
Coulomb interactions. 



P{Xmax = 2.428), where the jump in the pair-correlation 
occurs at the same t' c , indicating that charge mobility is 
due to pair motion. The results of Fig. 2 indicate that 
coexistence of CO and SC occurs over a very narrow re- 
gion of the parameter space, if at all, for this moderately 
strongly correlated case. This is in contrast to the mixed 
CO-SC state that is obtained when instead of t' and V 
the carrier concentration n is varied^. 

Within Eq. 1, the amplitude of the CO at if = in- 
creases with both \U\ and V . As might be expected, the 
stronger the CO in the t' = limit, the larger is the t' c at 
which the CO-to-SC transition occurs. Furthermore, it 
is known that for V = the CO and SC states are degen- 
erate in the square lattice. It is then to be expected that 
the transition is second order for weak V , independent 
of the value of \U\. We show a t' — \U\ phase diagram 
for fixed V = 1 in Fig. 3(a). The transition remains first 
order like for this moderate V until \U\ is large. Our 
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FIG. 3: Phase diagrams of the frustrated EHM at zero tem- 
perature in (a) the t' — \U\ plane for fixed V = 1; (b) the 
t' — V plane for fixed U = —4. The transition is first order 
(continuous) for large (small) V in (b). The SC phase persists 
in the t' > 1 region until t' S> 1. 



phase diagram in the t' — V space in Fig. 3(b) is for fixed 
1 17 1 = 4. Surprisingly, when both \U\ and V are large 
(|J7| = 4, V = 2) the ground state continues to be the 
checkerboard CO for t' as large as 0.9. 

The SC phase within Eq. Q]is robust and occurs over a 
broad region in the phase diagrams in Figs. 3(a) and (b). 
For the \U\ values of interest SC is lost only for t' 3> 1. 
This is in sharp contrast to the results obtained within 
the spin frustration model with repulsive U, where even 
within the approximate methods that claim SC, the bulk 
of the t' — U parameter space is occupied by AFM and 
PM phases, and SC occurs as an intermediate phase for 
a narrow range of t' for each U-. 

We now speculate how the negative-U EHM in 
the weak \U\ limit may be relevant to the CTS, (3- 
Nao.33V205 and other n = \ systems. Our fundamental 
premise is that the behavior of EtMe3Z[Pd(dmit)2]2 and 
k— ET 2 Cu2(CN) 3 with nearly isotropic lattices is repre- 
sentative of the more anisotropic systems at very low 
temperatures and under pressure, when V' and t' are 
large. The mapping of the molecular n = 5 lattice into 
the effective n = 1 repulsive Hubbard Hamiltonian that 
describes the AFM phase requires homogeneous charge 
population on the dimer unit cells^. We postulate that 
under pressure the AFM phase switches over to a static or 
fluctuating VBS with inhomogenous charge distribution 
(see below), and the effective n = 1 Hamiltonian that de- 
scribes this state is different. We rationalize this hypoth- 
esis based on our recent work. The complete Hamiltonian 
for these systems must start from the n = i repulsive-U 
EHM. In a series of papers^, we have established that the 
ground state of the n = 5 repulsive Hamiltonian is often 




FIG. 4: (a) BCDW corresponding to the VBS in the 2D or- 
ganic layer. Filled (unfilled) ellipses correspond to charge-rich 
(charge-poor) molecules, (b) Pairs of charge-rich and charge- 
poor molecules constitute the double occupancy (filled circle) 
and vacancy (empty circle) in the effective n — 1 lattice. 



a Bond-Charge-Density Wave (BCDW) , with charge oc- 
cupancy • ■ • 1100 • • • , where T' and '0' refer to molecular 
charges 0.5 + e and 0.5 — e, respectively. This is a quan- 
tum effect driven by the antiferromagnetic correlations 
due to the repulsive U and dominates over the classical 
effect due to V that favors the formation of the Wigner 
crystal • • • 1010 • • • for V < V C (U), where V C (U) can be as 
large as 3 for realistic U—. The BCDW is enhanced by 
the intersite electron-phonon interactions that modulate 
the hopping integral, and the intrasite Holstein electron- 
molecular vibration couplings^ -. Most importantly, even 
when at high temperatures there occurs a Wigner crys- 
tal CO in the (TMTTF) 2 X, the low temperature spin- 
Peierls phase is the BCDW—. We postulate that the 
• • • 1100 • • • is a bipolaron density wave that can be mod- 
eled by an effective n = 1 negative-U EHM, where the 
effective sites are alternating pairs of occupied (1-1) and 
unoccupied (0-0) molecules, respectively. That the sites 
of a negative-U EHM for a complex system can be com- 
posite ones has been suggested previously^. Note that 
the traditional mapping of n = 5 into effective n = 1 
AFM 7 also assumes composite dimer sites, with the only 
difference that our proposed mapping is valid for the in- 
homogenous charge distribution. 

The ■••1100--- charge occupancy persists for di- 
mensionality greater than 1—. Experimentally, 
the 2D BCDW has been seen in the weakly 2D 
(TMTSF) 2 X, the triangular lattice 6»-ET 2 X, the lad- 
der materials (DTTTF) 2 M(mnt) 2 (M=Au, Cu), and 
EtMe3[Pd(dmit) 2 ]2- The co-existing 2kp charge- and 
spin-density waves^ in (TMTSF^X has been explained 
as a Bond-Charge-Spin Density Wave^. The exper- 
imentally observed horizontal charge stripe in the a- 
ET 2 X and 6»-ET 2 X consists of •••1100- •• charge oc- 
cupancies occurring along the two directions of strong 
hoppingsi^. The large spin gap in the two-chain 
(DTTTF) 2 M(mnt)2^- can be understood within the 
n = h zigzag ladder modelSS, within which there oc- 
cur interchain singlets and the • • • 1100 • • ■ charge occu- 
pancy along the zigzag diagonals^. We point out that 
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charge occupancies and the bonding patterns in the VBS 
EtMe 3 P[Pd(dmit) 2 ]2 (see Fig. 3(b) in Ref. li and Fig. 4) 
correspond precisely to the ■ ■ ■ 1100 • • ■ BCDW. Finally, 
/3-Nao.33V205 consists of n = | V chains and ladders, 
and there is extensive literature on n.n. intersite bipo- 
larons in this and related systems^ (the earlier literature 
did not emphasize the n = \ carrier concentration that 
is crucial for the mapping into the effective n = 1 Hamil- 
tonian) . 

We therefore postulate that n = | systems at low tem- 
peratures and under pressure can be understood qual- 
itatively within the effective n = 1 negative-?/ EHM. 
This "mapping" is seen in Fig. 4, where we have shown 
the intersite bipolaron CDW corresponding to the VBS 
in EtMe3P[Pd(dmit)2]2*^, as well as the corresponding 
n = 1 checkerboard CO. We speculate that for small t' 
within the n = i molecular system, the occupancies of 
the dimer unit cells are indeed homogeneous, which gives 
the AFM in this region. The larger t' or lattice distor- 
tion (or both) at low temperatures give the insulating 
VBS with CO, and at still larger t' the intersite spin sin- 
glets become mobile, which is approximately modeled in 
our calculations within the effective negative-?/ EHM. We 
emphasize that intersite (as opposed to intrasite) bipo- 
larons have been shown to be light and mobile, especially 
in the triangular lattice^, and will be even more so here 



as the binding is due to antiferromagnetic correlations 
and not overscreening of V by electron-phonon interac- 
tions. Finally, although SC within the effective Hamil- 
tonian is s-wavc, this need not be true within the actual 
n = o repulsive-C/ EHM with intersite pairs. 

In summary, we have shown that it is indeed possi- 
ble to have frustration-driven SC, within a model that 
emphasizes charge as opposed to spin-correlations. We 
have also suggested that the n = \ VBS insulator can be 
conceived as an effective n — 1 bipolaron CDW. A com- 
plete theory of SC in the CTS should of course demon- 
strate the CO-to-SC transition within the actual n = i 
repulsive-?/ EHM with electron-phonon interactions, but 
we believe that our work on the effective n = 1 Hamil- 
tonian gives the insight necessary for constructing such 
a theory. In particular, the approach opens up the pos- 
sibility of describing all CTS as well as inorganic n = | 
systems within one simple model. An interesting aspect 
of our work is that the configuration space pairing in the 
VBS, and by implication in the superconducting state 
within the present model occurs due to a co-operative in- 
teraction between AFM correlations and electron-phonon 
interactions, and thus SC mediated by these two interac- 
tions need not be mutually exclusive. 

This work was supported by the Department of Energy 
grant DE-FG02-06ER46315. 
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